J/ip Production vs Centrality, Transverse Momentum, and Rapidity 
in Au+Au collisions at ^/snn = 200 GeV 
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The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured J /if) pro- 
duction for rapidities —2.2 < y < 2.2 in Au + Au collisions at ^/snn = 200 GeV. The J /if) invariant 
yield and nuclear modification factor Raa as a function of centrality, transverse momentum and 
rapidity are reported. A suppression of J /if) relative to binary collision scaling of proton- proton 
reaction yields is observed. Models which describe the lower energy J/ip data at the Super Proton 
Synchrotron (SPS) invoking only J /if) destruction based on the local medium density would predict 
a significantly larger suppression at RHIC and more suppression at mid rapidity than at forward 
rapidity. Both trends are contradicted by our data. 

PACS numbers: 25.75.Dw 



The Quark-Gluon-Plasma (QGP) is a state of decon- 
fined quarks and gluons which is predicted by lattice 
Quantum Chromodynamics (QCD) calculations to be 
formed above a temperature of order T c = 170 MeV for 
a baryon chemical potential /if, = [l[ . Heavy quarkonia 
(J/ip, if)', Xc and T) have long been considered one of 
the most promising probes to study formation and prop- 
erties of QGP. In the deconfmed state, the attraction 
between heavy quarks and anti-quarks is predicted to be 
reduced due to dynamic screening effects, leading to the 
suppression of heavy quarkonia yield. The strength of 
the suppression depends on the binding energies of the 
quarkonia and the temperature of the surrounding sys- 
tem [2|. Recent lattice QCD calculations suggest that 
the J/ip may not dissociate until well above T c [1, 0|- 
On the other hand \c and if)' which contribute to the to- 
tal J/ip yield via decay are expected to dissolve at lower 
temperatures due to smaller binding energies. 

A J/ip suppression observed at lower energies by the 
NA50 experiment at the SPS 0, @ could be reproduced 
by various theoretical calculations, some invoking QGP 
formation and others not. A larger suppression is ex- 
pected at RHIC compared to SPS due to the larger en- 
ergy density of the medium created @, HI • On the other- 
hand, several models predict that the J/ip yield will re- 
sult from a balance between destruction due to thermal 
gluons and enhancement due to coalescence of uncorre- 
cted cc pairs 0, , which are produced abundantly in 
the initial collisions at RHIC energy [ill ]. Cold nu- 
clear matter (CNM) effects such as nuclear absorption, 
shadowing and anti-shadowing are also expected to mod- 
ify the J/ip yield. PHENIX d + Au data show that CNM 
effects are smaller at RHIC than those observed at lower 
energy [lji] and can be reproduced by a nuclear absorp- 
tion cross-section of up to 3 mb plus nuclear shadow- 
ing M- 

We report results on J/ip production measured by the 
PHENIX collaboration at mid-rapidity (\y\ < 0.35) via 



e + e~ decay and at forward rapidity (\y\ £ [1.2,2.2]) via 
fi + fi" decay in Au + Au collisions at ^/snn = 200 GeV. 
These results do not separate primordial J/ip and J/ip 
from \c, ip' or B decay. The J/ip invariant yields as 
a function of centrality, rapidity (y) and transverse mo- 
mentum (pt) are shown. They are combined with the 
J/ip yield measured in p + p collisions [3] to form the 
J/ip nuclear modification factor Raa- 

The PHENIX apparatus is described in [HI]. At mid- 
rapidity electrons are measured with two spectrometers 
consisting of Drift Chambers (DC), Pad Chambers (PC), 
Ring Imaging Cerenkov Counters (RICH), and Electro- 
magnetic Calorimeters (EMCal). They are identified by 
matching charged tracks reconstructed with the DC and 
the first layer of the PC to clusters in the EMCal and to 
hits in the RICH. The energy-momentum matching re- 
quirement is (E/p — 1) > —2.5 standard deviations (a). 
The position matching between the track and the clus- 
ter in the EMCal is < 2.5<r (4(7 ) in azimuth and along 
the beam axis, for central (peripheral) collisions. For the 
RICH, at least 4 (2) matching hits are required. Muons 
are measured with two forward spectrometers consist- 
ing of a front absorber to stop most hadrons produced 
in the collision, cathode strip chambers (MuTr) which 
provide momentum information and a Muon Identifier 
(MuID) which uses alternating layers of steel absorber 
and Iarocci tubes. Charged particle trajectories are first 
reconstructed in the MuID then in the MuTr. They must 
reach the last plane of the MuID and have a good geo- 
metrical match between the MuID and the MuTr to be 
identified as muons. The matching is < 9° for the slope 
and < 15 (20) cm for the position in the first layer of 
the MuID at positive (negative) rapidity. The collision 
centrality is determined using two Beam-Beam Counters 
(BBC) and Zero Degree Calorimeters (ZDC) 0. 

The data used for this analysis were collected during 
the 2004 run at RHIC using a minimum bias trigger (a 
coincidence of the two BBC) which covers 92 ±3% of the 
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Au + Au inelastic cross-section. After quality assurance 
and vertex cut {\z\ < 30 cm), 9.9 x 10 s (1.1 x 10 9 ) events 
were analyzed for mid (forward) rapidity, corresponding 
to an integrated luminosity of 157 fib -1 (174 /ib -1 ). The 
forward rapidity data were filtered using an offline lcvel-2 
trigger which provides a fast reconstruction of the parti- 
cle trajectory in the MuID. Events are accepted by this 
filter when at least two good quality tracks reaching the 
last plane of the MuID are found within the acceptance. 

The J/tp yield is obtained from the unlike sign dilep- 
ton invariant mass distribution after subtracting the 
combinatorial background using an event-mixing tech- 
nique. The background is normalized to the real data 
using the like-sign dilepton invariant mass distribution, 
2VN++N — , with N++ (N — ) being the number of 
positive (negative) dilepton pairs. The accuracy of the 
normalization is estimated to be 2 % and accounted for 
in the systematic errors. At midrapidity the J/tp mass 
resolution is ~ 35 MeV/c 2 . The number of J/tp is deter- 
mined by counting the remaining unlike sign pairs in the 
mass range 2.9 < M e + e - < 3.3 GeV/c 2 . This number is 
corrected by the estimated contribution of the dielectron 
continuum and the loss due to the radiative tail. A total 
of ~ 1000 J/tp are obtained and the signal to background 
(S/B) varies from 0.5 for central collisions to 15 for pe- 
ripheral collisions. At forward rapidity, the J/tp mass res- 
olution varies from 150 to 200 MeV/c 2 and is larger than 
at midrapidity primarily because of the multiple scatter- 
ing and energy loss straggling in the front absorber. The 
residual background (notably from the open charm pairs 
and Drell-Yan) in the unlike-sign invariant mass distri- 
bution is evaluated using an exponential form. The J/tp 
signal is estimated with direct counting of the remaining 
pairs in the mass range 2.6 < M^^- < 3.6 GeV/c 2 and 
using a fit with different line shapes. The average of the 
resulting values is used as the number of J/tp and their 
dispersion is included in the systematic error. A total 
of ~ 4500 J/tp are obtained and S/B varies from 0.2 for 
central collisions to 3 for peripheral collisions. 

The J/tp invariant yield in a given centrality, p^ and 
y bin is: 

Bg d 2 N J/lp = 1 N JH , 
2irpT dp^dy 2irpT N evt AyApTAs 

with Bu being the branching ratio for J/tp — > Nj/^ 
the number of J/tp measured in the centrality bin; N ev t 
the corresponding number of events; Ay the rapidity 
range; Apt the transverse momentum range and Ae the 
acceptance and efficiency correction for J/tp. Ae is deter- 
mined by full GEANT simulation. It decreases with the 
collision centrality due to overlapping hits in the RICH, 
EMCal and MuTr, leading to an increasing amount of 
mis-reconstructed tracks which are then rejected by the 
analysis cuts. This effect is evaluated by embedding sim- 
ulated J/tp in real events. For the most central collisions 



TABLE I: Sources of systematic errors on the J/tp invariant 
yield. Columns 2 (3) are the average values at mid (forward) 
rapidity. When two values are given, the first (second) is for 
peripheral (central) collisions. Errors of type A (type B) are 
point to point uncorrelated (correlated). 



source 


\y\ < 0.35 


\y\ € [1.2,2.2] 


type 


signal extraction 


6.5 to 9 % 


4 to 24 % 


A 


acceptance 


6 % 


10 % 


B 


efficiency 


4.5 to 8 % 


4 to 16 % 


B 


run by run variation 


4 % 


5 % 


B 


input y, pt distributions 


2 % 


4 % 


B 



the efficiency loss is 20 % at mid rapidity and 75 % (50 %) 
at positive (negative) rapidity. 

The nuclear modification factor in a given centrality, 
Pt and y bin is: 

d 2 Nj^/dp T dy 
UAA N coll d 2 N p / /ip /dp T dy (Z) 

with d 2 NjA/dpTdy being the J/tp yield in Au + Au col- 
lisions in the centrality bin, N co \\ the corresponding mean 
number of binary collisions and d 2 N^j^/ dp^dy the J/tp 
yield in p + p inelastic collisions. 

The systematic errors on the J/tp invariant yield (Ta- 
ble H| are grouped into three categories: point to point 
uncorrelated (type A) for which the points can move in- 
dependently one from the other; point to point correlated 
(type B) for which the points can move coherently though 
not necessarily by the same amount; and global errors for 
which all points move by the same relative amount. Sta- 
tistical and uncorrelated systematic errors are summed 
in quadrature and represented with vertical bars; corre- 
lated systematic errors are represented with boxes and 
different colors/symbols are used for forward and mid 
rapidity because they are independent; global systematic 
errors are quoted directly on the figures. For Raa, addi- 
tional systematic errors are associated with uncertainties 
in the calculation of iV C oii (10 to 28 %) and the J/tp yield 
in p + p (12 % and 7 % at mid and forward rapidity, 
respectively). On the other hand the systematic errors 
that are common to Au + Au and p + p cancel. 

Figure Q] shows the J/tp yield vs. px for different cen- 
trality bins (see Table ITT1 for the corresponding number of 
participants, N part ). Data from the two muon spectrom- 
eters are combined to obtain the forward rapidity points. 
In each centrality bin, the J/tp mean square transverse 
momentum, ^px 2 ), is numerically calculated from the 
data for pt < 5 GeV/c and is shown in Table HT1 The 
first error corresponds to the statistical and uncorrelated 
systematic error on the J/tp yield. The second corre- 
sponds to the correlated systematic error. At midrapid- 
ity the (pt 2 ) shows no variation versus centrality within 
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FIG. 1: J/tp invariant yield vs. pt for different centrality 
bins in Au ± Au collisions and in p + p collisions 14] . The 
left (right) panel corresponds to mid (forward) rapidity. See 
text for description of the errors. 



FIG. 2: J/tp invariant yield vs. y for different centrality bins 
in Au ± Au collisions and for p + p collisions. Open (filled) 
circles are for mid (forward) rapidity Au + Au data. Black 
squares are for p + p data [1J|. See text for description of the 
errors. 



TABLE II: Characterization of the J/ip pt and y distribu- 
tions. Column 3 (4): J/tp (pt 2 ) calculated for Pt < 5 GeV/c 
at mid (forward) rapidity for different centrality bins in 
Au + Au collisions and in p + p collisions. The errors are 
descibed in the text. Column 5: Calculated RMS of the cor- 
responding J/tp y distributions. 

cent iV part (pt 2 ) (GeV/c) 2 (p T 2 ) (GeV/c) 2 y RMS 

(%) ly|<0.35 1.2 < |y| < 2.2 

0-20 280 3.6 ±0.6 ±0.1 4.4 ±0.4 ±0.4 1.32 ±0.06 

20-40 140 4.6 ±0.5 ±0.1 4.6 ± 0.3 ± 0.4 1.30 ±0.05 

40-60 60 4.5 ±0.7 ±0.2 3.7 ±0.2 ±0.3 1.40 ±0.04 

60-92 14 3.6 ±0.9 ±0.2 3.3 ± 0.3 ± 0.2 1.43 ±0.04 

p + p 2 4.1 ±0.2 ±0.1 3.4±0.1±0.1 1.41 ± 0.03 



the error bars. It increases slightly with -/V part at forward 
rapidity. 

Figure [2] shows the J/tp yield (integrated over pr) vs. 
y for different centrality bins. The root mean square 
(RMS) of each distribution is shown in Table HT1 For the 
two most peripheral bins (40-60 % and 60-93 %) the RMS 
is compatible with that measured in p + p collisions. For 
the most central bins (0-20 % and 20-40 %), the RMS is 
smaller by about 2a. 

Figures [3] and [4] show the J/tp Raa vs. pt and y, re- 
spectively, for different centrality bins. Figure[5^a) shows 
the pt integrated Raa vs. 7V par t at mid and forward ra- 
pidity. For each rapidity, Raa decreases with increasing 



For the most central collisions, Raa is below 0.3 
(0.2) at mid (forward) rapidity. Figure G^b) shows the 



part ■ 



ratio of forward/mid rapidity Raa vs. A^ part . The ra- 
tio first decreases then reaches a plateau of about 0.6 for 

iV part > 100. 

We observed a significant J/tp suppression relative to 
binary scaling of proton-proton is observed for central 
Au ± Au collisions at RHIC. The magnitude of the sup- 
pression is similar to that observed at the SPS Q and 
greater than the suppression expected by extrapolating 
the cold nuclear matter effects measured in d + Au colli- 
sions [lH, [l3j]. Models that describe the SPS data using 
a J/tp and/or \c an( i ty' suppression based on the local 
density predict a significantly larger suppression at RHIC 
than SPS and more suppression at mid rapidity than 
at forward rapidity 0, Q . Both trends are contradicted 
by our data. Additionally, the J/tp mean square trans- 
verse momentum, restricted to px < 5 GeV/c, shows 
little dependence on centrality. Various models of J/tp 
production and suppression, which predict very different 
transverse momentum and rapidity dependencies, can be 
significantly constrained by the data presented here and 
recent results on open charm [il| . 
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FIG. 3: J/ip Raa vs. pt for several centrality bins in Au + Au 
collisions. Mid (forward) rapidity data are shown with open 
(filled) circles. See text for description of the errors. 
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FIG. 4: J/ip Raa vs. y for different centrality bins. Open 
(filled) circles are for mid (forward) rapidity. See text for 
description of the errors. 
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FIG. 5: (a) J/ip Raa vs. Apart for Au + Au collisions. Mid 
(forward) rapidity data are shown with open (filled) circles, 
(b) Ratio of forward/mid rapidity J/ip Raa vs. Af par t. For 
the two most central bins, mid rapidity points have been com- 
bined to form the ratio with the forward rapidity points. See 
text for description of the errors. 



